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COMPUTATION  OF  THE  NONEQUILIBRIUM  FLOW  IN  A  GASDYNAMIC  LASER 

* 

Yen  Hai-hsing,  Chen  Li-yin 

(Institute  of  Mechanics,  Academia  Sinica) 

ABSTRACT 

In  this  paper,  we  present  a  three-modes-four- 
temperatures  vibrational  relaxation  model  for  the  C02~ 
N2~H2°  laser  system,  and  give  a  set  of  fairly  rigorous 
relaxation  equations .  We  have  analyzed  a  series  of 
problems  related  to  computations  of  pseudo-one-dimen¬ 
sional  non-equilibrium  flows,  and  have  carried  out 
massive  numerical  computations  using  the  data  we  obtained 
for  the  rate  of  relaxation.  Our  computed  results,  unlike 
those  of  other  authors,  agree  fairly  well  with  experi¬ 
mental  results.  The  performance  of  the  laser  can  be 
effectively  improved  by  increasing  the  expansion  area 
ratio.  Appropriately  decreasing  the  height  of  the  throat 
also  proves  to  be  helpful.  For  combustion  type  gas- 
dynamic  lasers,  there  is  an  optimal  value  for  the  stag¬ 
nation  temperature,  which  lies  between  1400  and  1600 °K. 
Contrary  to  Anderson's  conclusion,  our  computed  results 
show  that  in  a  nozzle  flow  with  a  high  stagnation  temper¬ 
ature  and  a  large  area  ratio,  the  optimal  water  content  is 
still  in  the  neighborhood  of  1% .  The  performance  of  the 
unit  quickly  deteriorates  with  increasing  water  content. 

We  have  studied  the  effect  of  the  various  section  of  the 
shape-curve  of  the  nozzle  on  the  performance  of  the  unit. 
We  have  also  given  a  preliminary  discussion  of  the  effect 
of  the  relaxation  model,  equations  and  data  on  the  compu¬ 
tation,  and  have  shown  that  the  results  of  computation  are 
very  much  affected  by  these  factors. 


I .  Introduction 


Gasdynamic  lasers  are  a  new  type  of  high  power  laser  that 
appeared  near  the  end  of  the  60' s.  These  lasers  have  received  great 
attention  in  recent  years,  and  have  developed  very  rapidly  [1] . 
Theoretical  analyses  have  played  an  important  role  in  every  stage 
of  the  development  of  gasdynamic  lasers.  In  particular,  computa¬ 
tions  of  pseudo-one-dimensional  vibrational  non-equilibrium  flow 
have  provided  useful  means  for  designing  and  improving  gasdynamic 
lasers . 

The  computation  of  pseudo-one-dimensional  non-equilibrium  flow 
falls  into  two  types.  One  is  the  steady-state  method  first  adopted 
by  ,  et  al.,  [23  —  [6]  .  In  this  approach,  the  pseudo-one- 

dimensional  flow  equations  are  solved  simultaneously  with  the  vibra¬ 
tional  relaxation  equations.  The  set  of  ordinary  differential  equa¬ 
tions  obtained  are  solved  by  means  of  the  Runge-Kutta  method.  As 
there  are  singularities  at  the  sonic  points  near  the  convergent- 
divergent  nozzle,  it  is  not  possible  to  find  a  simple  unified  solu¬ 
tion  for  the  entire  nozzle.  To  simplify  the  calculations,  the  sub¬ 
sonic  section  of  the  nozzle  is  usually  assumed  to  be  in  equilibrium, 
and  the  equilibrium  flow  is  calculated  starting  from  the  throat.  To 
overcome  the  difficulty  with  the  singularities,  certain  techniques 
of  parameter  adjustment  are  used  to  obtain  convergent  solutions  for 
the  set  of  equations. 

The  other  type  is  the  non-steady-state  method  adopted  by  Anderson 
[7],  In  this  approach,  a  distribution  of  initial  values  at  time  t=0 
is  first  assumed  for  each  of  the  unknown  variables  in  the  set  of  non¬ 
steady-state  equations.  Then,  a  difference  scheme  is  used  to  find 
the  values  of  the  variables  at  time  t+At  (the  time  step  size  At 
being  required  to  meet  certain  stability  conditions)  from  their  values 
at  time  t.  The  process  is  repeated  until,  at  a  fairly  large  t,  all 
the  parameters  converge  to  their  steady-state  values .  These  types  of 
equations  are  hyperbolic  for  the  entire  nozzle.  No  singularity 


exists  at  the  throat  region.  The  entire  nozzle  can  thus  be  solved 
for  as  a  whole.  As  a  result,  the  nonequilibrium  effects  at  the  sub¬ 
sonic  section  of  the  nozzle  are  automatically  included  in  the  com¬ 
putation,  and  possible  errors  due  to  artificially  introduced  adjust¬ 
ments  in  the  numerical  calculations  are  avoided.  The  obvious  advan¬ 
tage  of  this  method  has  made  it  very  popular  in  theoretical  computa¬ 
tions  on  gasdynamic  lasers  [8]  —  [13]  . 

However,  researchers  abroad  have  either  used  simplified  relax¬ 
ation  models  in  their  computations,  or  employed  erroneous  relaxation 
equations.  In  addition,  most  of  these  computations  have  not  included 
the  latest  relaxation  data.  Hence,  certain  deviations  or  contradic¬ 
tions  exist  between  their  computed  results  and  the  experimentally 
measured  values  [14] , [10] . 


In  this  paper,  we  present  a  three-modes-f our- temperature  relax¬ 
ation  model,  and  give  a  set  of  fairly  rigorous  relaxation  equations. 
The  equations  for  computing  nonequilibrium  flows  will  be  listed, 
some  problems  related  to  the  numerical  calculations  will  be  discussed 
and  the  main  computed  results  will  be  given  and  compared  with  the 
theoretical  and  experimental  results  obtained  by  others .  We  will 
also  give  an  analysis  of  certain  related  problems. 

II .  Relaxation  model  and  relaxation  equations 

The  rate  of  vibrational  relaxation  of  the  "catalyst"  H20  in  the 
CO2“N2-H20  laser  system  is  very  high,  so  that  the  various  vibrational 
modes  can  always  be  regarded  as  being  approximately  in  equilibrium 
with  the  average  motion.  Thus,  the  V-V  exchange  among  the  molecules 
of  H20  and  those  of  C02  and  N2  is  equivalent  to  the  result  of  regard¬ 
ing  H20  as  collision  partners  for  C02  and  N2  molecules  in  their  mole¬ 
cular  relaxation  process  [15]  .  The  vibrational  energy  level  diagram 
and  energy  exchange  processes  for  the  laser  system  are  given  in 
Figure  1.  The  vibrational  relaxation  processes  are: 


V-T  exchange 


intermo lecular 
V-V  exchange 

intramolecular 
V-V  exchange 
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In  the  .above,  *  denotes  a  quantum  of  vibrational  energy.  The  colli¬ 
sion  partner  M  —  CO,,  N„  HjO(C,N,H).  .  The  relaxation  equations  are 
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In  the  above,  E^,  E2,  ,  EN  are, 

stored  in  the  vibrational  modes 
unit  mass  (g)  of  the  gas  mixtures 


respectively  the  vibrational  energy 
COj(vi) ,  COjQvi),  CO,(vJ  and  N2  of  a 
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Here  =  N2  =  N3  =  Nc  and  NN  are,  respectively,  the  number  of  mole¬ 
cules  of  C02  and  N2  per  gram  of  the  gas  mixture:  g.  is  the  degree  of 
degeneracy,  which  is  1  except  for  g2  which  is  2;  is  the  vibra¬ 
tional  temperature;  9^  is  the  characteristic  vibrational  temperature; 
k  is  Boltzmann's  constant. 


Figure  1.  C02-N2  vibrational  energy  level  diagram 

1 — vibrational  model  I;  2 — vibrational  mode  II;  3 — vibrational  mode  III 


It  can  be  derived  from  the  general  vibrational  relaxation  equations 
[16]  that 
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In  the  above,  the  equilibrium  vibrational  energy  E2(T)  and  E^ (T)  have, 
respectively,  been  obtained  from  equation  (10)  by  substituting  the 
vibrational  temperature  T-  and  T  with  the  temperature  T  of  the  aver- 

/«r\  W 


age  motion .  N 


is  the  density  of  particles  of  the  constituent  M 


per  unit  volume  (cm-3):  NIT  —  #*V„ 

(16) 

p  is  density  (g/cm^) ,  and  K  is  the  coefficient  of  the  rate  of  the 
vibrational  relaxation  process . 

As  Fermi  resonance  exists  between  the  v«»  modes  of  C02,  the 
process  of  CO, (*,-+  2»a)  is  very  rapid.  In  other  words,  the  value  of 
k[Cl(l -* 0)7  C2(0 -* 2)  is  very  large.  If  at  a  certain  instant  the 
vibrational  energy  of  the  >4  and  v,  vibrational  modes  does  not  satisfy 

(-&- + 0  M«-Y  -  (-M  (-5t_  +  iY  *  c  n7, 

\NcHd ,  AlVcW  V.Ve^,/  \2.Vc^,  )  d7) 

then  it  can  be  seen  from  equation  (15)  that  a  large  value  of  K  will 
result  in  a  high  rate  of  change  with  time  of  E^ ,  and  the  *» 

vibrational  mode  will  equilibrate  with  the  h  mode  in  a  very  short 


time.  This  means  that  as  long  as  and  do  not  satisfy  equation 
(17),  the  rapid  CO,(*, -* relaxation  process  will  cause  the  left- 
hand  side  of  equation  (17)  to  go  to  zero  in  a  negligible  interval 
of  time.  Thus,  we  think  equation  (17)  is  a  very  good  approximation. 
It  can  be  easily  derived  from  equation  (17)  that 

4  -  igl  —  gi_=Jg?  (18) 

T,  T ,  T 

This  is  the  condition  for  equilibrium  of  the  *’>>  »»  modes.  It  is 
worth  noting  that  the  condition  T^  =  T2  used  by  Anderson,  et  al., 
is  valid  only  when  one  further  assumes  that  0,  «=»  20,  or  =  T.  In 
all  other  cases,  >  T^  >  T. 


Now  we  can  combine  the  modes  to  form  the  vibrational 

mode  I,  with  vibrational  energy  given  by 

E'-E.  +  E.  (19) 

We  have  from  equations  (6)  and  (7) 

,  0,-2 0,/iEA  +  (ih\  +  (4Ei\ 

dt  \<t*  0,  \  it  /»,-!», +  V  it  /»,-**.  w*  •  ».-T  (2°) 

radiative 

Let  us  analyze  the  magnitude  of  each  term.  Although  the  terms 
enclosed  in  the  right-most  sauare  brackets  in  equation  (15)  are  very 
small,  as  K  is  fairly  large,  is  not  necessarily  small.  In 

fact,  in  the  absence  of  a  radiation  field,  we  have 


("0 

(«l\  +  ~  C  +  1  )* 

'it  ) »i-i  V dt  1 


with  values  equal  to  0.-0,  0.39  and  0.45,  respectively,  when  T2  is 

350,  1000  and  1500°K.  In  other  words,  it  is  not  appropriate  to 

eliminate  the  term  from  equations  (6)  and  (7)  .  However, 

'it  /»,->•  g.  —  7ejdE,\ 

in  equation  (20),  as  (0,  —  20, )/0,  —  0.0387  ,  the  term  »* 

is  negligible  compared  to  the  other  terms  in  this  equa¬ 
tion.  A  confession  is  due  here  in  that  part  of  the  reason  for  want¬ 
ing  to  eliminate  this  term  lies  in  the  lack  of  reliable  data  for 
K(Cl(l-*0),  C2(0 -* 2)  ;  M],  Otherwise,  one  could  use  equation  (15) 

directly  in  equations  (6)  and  (7),  and  compute,  without  much  diffi¬ 
culty,  according  to  a  four-vibrational  modes-four-temperatures  model 


We  can  use  the  equilibrium  condition  equation  (18)  to  transform 
equation  (20)  into 


This  equation,  along  with  equations  (8)  and  (9) ,  make  up  the  set  of 
equations  for  the  three-modes-four-temperatures  relaxation  model  of 
this  paper.  The  various  terms  in  the  above  equations  are  given  by 
equations  (11)-(14). 

The  inaccuracy  of  the  data  for  relaxation  rate  has  a  large 
effect  on  the  outcome  of  nonequilibrium  flow  computations  [17,13]. 
Most  authors  have  based  their  computations  on  the  data  given  in  [15]. 
Many  experimental  results  on  relaxation  were  later  published,  and 
these  contained  additional  data  or  modifications  of  previous  data 
that  have  a  rather  great  effect  on  particle  number  inversion.  We 
have  collected  and  processed  experimental  data  of  the  relaxation 
rate  coefficients  for  the  processes  given  in  equations  (l)-(4),  and 
have  arrived  at  the  following  relation  for  K  and  temperature  T: 

(“ ■)  —  A  +  BT ~xn  +  CT~tn  +  DT~' 

The  coefficients  A,  B,  C  and  D  have  been  given  for  each  of  the  reac¬ 
tions.  G  is  a  combinatorial  factor.  For  reactions  ( 1) — ( 3) ,  G  =  1; 
for  reaction  (4),  G  =  4. 

The  computations  made  in  this  work  are  based  on  the  above  three- 
modes-four-temperatures  model  and  the  more  rigorous  set  of  relaxation 
equations,  and  make  use  of  the  relaxation  data  that  we  hhve  collected 
and  processed  ourselves. 

The  relaxation  model  and  equations  of  Anderson  are  still  being 
used  widely.  A  brief  discussion  of  these  is  in  order  here.  Ander¬ 
son  regards  CO,(*>,)  to  be  in  equilibrium  with  the  vibrational  modes 


of  N2/  and  considers  their  vibrational  temperatures  to  be  the  same, 
thus  forming  mode  II.  This  kind  of  assumption  is  not  permissible 
in  the  presence  of  radiation  fields .  As  the  lasing  process  depends 
to  a  large  degree  on  the  relaxation  rate  of  the  transfer  of  the 
vibrational  energy  stored  in  N2  to  the  c0^V!-)  vibrational  mode, 
regarding  the  two  vibrational  modes  as  being  in  equilibrium  is  equi¬ 
valent  to  considering  the  relaxation  rate  to  be  infinitely  large. 
This  is  obviously  not  appropriate.  Even  in  calculations  involving 
small  signal  gains,  errors  will  result.  For  the  typical  laser  com¬ 
position  XC:XS’XU  *■  0.14:0.85:0.01  ,  the  computed  relaxation  times  are 
given  in  Table  1. 


TABLE  1 


T(*K) 

300  j 

i  * 

1000 

1500 

3.02 

2.03 

0.405 

0.096 

rrco,i 

0.078 

0.174 

0.480 

0.630 

It  can  be  seen  from  the  above  table  that  at  high  temperatures,  the 
relaxation  rate  of  COj(v0-*.V,  is  even  slower  than  that  of  COj(*j-*  3*0. 
Furthermore,  the  results  of  computation  show  that  for  a  nozzle  with 
A/Am*»  30,  and  h*  =  0.3mm,  T.- l500°K,ft- 30ann  XC:X,,:XH  -  0.14:0.85:0.01 . 

T^  is  higher  than  T3  by  up  to  130 °K  inside  the  nozzle;  and  at  the 
exit,  T^  is  also  higher  than  T^,  by  40°K.  This  difference  between 
TN  and  T3  is  even  more  pronounced  when  A/A*  is  increased  and  h*  is 
reduced . 


Another  problem  lies  in  the  fact  that  Anderson  used  the  linear¬ 
ized  relaxation  equations  for. both  modes  I  and  II: 

*£- —IE,(T) +  E,(T)- E,- E,1  (23) 

«  r».-7 

~-^(£,(r)  +  E„(r)-£,-E.v]  (24) 

Using  equation  (24)  is  equivalent  to  expanding  the  relaxation  term 
for  COj(*, -* 3vjj  ,  i.e.,  equation  (12),  and  taking  only  one  linear  term 
omitting  the  others.  This  is  only  valid  in  the  case  where  T2  =  T, 
which  is  not  true  in  gas  dynamic  lasers.  In  equation  (23),  Anderson 


has  omitted  the  effect  of  the  C0,(*>, -* 3v,)  process  on  mode  I,  and  has 
artificially  introduced  the  term  [£,(7)  —  E,]/rp><T  •  The  term 

E,(T)  —  e,  __  jg»  ± 

E,(T)  -  E,  (**•"•  +  l)(«vr  +  1) 

is  0.075  when  T2  %  T=300°K.  It  can  be  as  large  as  0.78  at  1000°K. 
Thus,  a  fairly  large  term  has  been  introduced  artificially.  Some 
of  Anderson's  computed  results  agree  better  with  experimental  results. 
This  may  be  due  to  the  fact  that  the  approximations  in  the  relaxation 
model  and  relaxation  equations  have  certain  errors  that  cancel  each 
other,  and  the  fact  that  the  not  very  accurate  relaxation  data  have 
compensated  the  inadequacies  of  the  model  and  the  equations. 

Ill .  Equations 


For  the  two-dimensional  nozzle  that  is  usually  used  in  gas- 
dynamic  lasers,  if  one  assumes  that  the  various  parameters  are  con¬ 
stant  across  the  flow,  the  problem  reduces  to  that  of  a  pseudo-one¬ 
dimensional  flow.  In  addition,  we  assume  the  fluid  to  be  an  ideal 
fluid,  and  neglect  any  inhomogeneity  in  the  flow  (such  as  boundary 
layer,  wake,  shock  waves,  etc.)  and  heat  transfer  effects.  Such 
approximations  work  better  for  a  single  nozzle  than  for  a  nozzle  array. 


After  p  is  eliminated  by  means  of  the  state  equation  p  =  pRT, 


the  pseudo-one-dimensional  flow 
mass  conservation  equation 

momentum  conservation  equation 

energy  conservation  equation 


equations  become: 
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dA 

—  0 

du 

(25) 

dt 

dx 
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dt 
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dx 

du  _ 

_  du  RT 

Sfi. 

-  R 

dT 

(26) 

dx 

P 

dx 

dx 

dE  _ 

•RT  dA 

RT 

du 

(27) 

dt 

A  dx 

dx 

In  the  above,  A  =  A(x)  denotes  the  cross-section  of  the  nozzle,  u 
is  the  flow  rate  and  E  is  the  total  internal  energy  of  1  g  of  gas 
mixture : 


2  2  tT\  —  1  « •  1 


R  is  the  gas  constant  for  1  g  of  gas  mixture,  is  the  gram-mole¬ 
cular  weight  of  H20,  and  t=4  denotes  the  N2  molecule. 

If  we  are  only  interested  in  computing  small  signal  gains,  then 
the  radiation  fields  are  considered  to  be  approximately  zero,  and 
one  has  ,  . 

V  it  , radia  ^  /radiative  mq\ 


Under  these  conditions,  we  can  write  equation  (27)  as 

dj  _  _u  37 _ 1 _ [3ft  -  6,/JEA 

dt  “  “  dx  (5  +  Xu)K  4-.  Wkdaal  ‘  3ft  \dt 

2  +  &  TV*"  -  IV 

+  (QA  |  +  +-EI^  +  RTp. 

'  dt  /»,«i  ft  '4*  /»,-»*  '  'nr*T  ^  0*  0* 

Equations  (25),  (26),  (30)  and  (9)  and  the  equation  (22)  obtained 
after  substituting  into  it  equation  (29)  make  up  a  complete  set  of 
equations,  with  p,  u,  T,  E2,  E^  and  E^  as  unknown  variables. 


After  the  steady-state  solutions  for  the  nonequilibrium  flow 
have  been  found  using  the  non-steady-state  method,  one  needs  to  fur¬ 
ther  compute  for  particle  number  inversion,  gain  and  the  maximal 
useable  laser  energy.  First,  we  find  the  vibrational  temperatures 

T1-T. .  Then,  we  obtain  the  density  of  particles  in  the  upper  and 
’  (y)  (w) 

lower  vibrational  energy  levels  of  C02,  N(qod  an<^  N(100)'  us^n9 

equatiohs  of  statistical  mechanics. 


The  gain  is  [18] 


C(.v)  -  ( NjV)  - 

v  *  ' 


.  .  ~  yCv) 

radiative 


In  the  above,  N, 


and  N. 


are,  respectively,  the  number  of  parti¬ 


cles  per  unit  volume  on  the  upper  and  lower  vibrational  energy  levels 
g^  denotes  the  degree  of  degeneracy,  is  the  central  frequency  of 
the  laser  transition,  c  is  the  speed  of  light,  g(v)  is  the  line 
shape  factor,  and  t  is  the  life  of  the  spontaneous  radiation. 


From  statistical  mechanics,  the  particle  number  inversion  of 
the  p(J)  transition  of  CO2 (001) -  (100)  is 


Zi 


-  Co-Vj  -  isrio  X  — ^ 


(32) 


with  the  factor 


0k»i  eip  I  +  DHL 


dMtxpl-6tMa  +  \)U  +  2)/T) 


(33) 


The  rotational  characteristic  temperature  is  determined  from  spectral 
constants  [19].  The  p(20)  transition  that  occurs  under  most  condi¬ 
tions  at  room  temperature  is 

«p(-235.5/T)  (34) 

0  —  0.992045exp(24.12/T)  (35^ 

When  T  is  1500,  1000,  500  and  300°K,  a  is  1.008,  1.016,  1.041  and 
1.075  respectively.  Obviously,  an  error  exists  in  the  formula 


&NZ2L-  -  (AT.V  -  AW)  «p(-234/T)  (36) 

that  Anderson  used  to  compute  gain.  When  N.jYl<<NnYl  ,  the  error  is 


(V) 


100  ”001' 


1-8%.  The  error  is  even  larger  when  N.  '  is  not  very  much  different 

(V)  i0U 

from  N00j_. 


Another  factor  in  the  gain  formula  equation  (31)  that  needs  to 
be  examined  is  the  line  shape  factor  g(v).  For  the  C02  laser  trans¬ 
ition,  the  natural  line  width  is  much  smaller  than  the  broadening  due 
to  collisions  and  the  Doppler  effect,  and  can  be  neglected.  At  high 
cavity  pressures,  the  line  width  is  mainly  due  to  collisions,  while 
at  low  cavity  pressures,  it  is  mainly  due  to  Doppler  broadening. 

For  a  laser  with  a  typical  composition  whose  cavity  is  at  room  tem¬ 
perature,  a  1%  error  will  result  from  considering  broadening  due  to 
collision  only,  when  the  cavity  pressure  is  at  0.1  atm,  while  a  3% 
error  will  result  if  the  same  is  done  when  the  cavity  pressure  is 
0.05  atm  [20].  When  only  broadening  due  to  collisions  is  taken  into 
consideration,  the  formula  for  calculating  the  gain  near  the  center 


of  the  spectral  line  is 


(37) 


<^00  - 

In  the  above  equation,  Avp  is  the  half-width  of  the  collision  broad¬ 
ening.  If  one  wishes  to  include  the  effect  of  Doppler  broadening, 
a  more  precise  formula  may  be  used  [20] . 

The  maximal  useable  laser  energy  defined  in  [7]  can  be  used  as 
a  rough  estimate  of  the  upper  limit  of  the  laser  energy  extractable 
from  the  lasing  medium.  It  is  the  difference  between  the  vibrational 
energy  stored  in  CC>2  (v^)  and  ^(v)  and  the  vibrational  energy  that 
still  remains  when  the  gain  has  dropped  to  zero,  multiplied  by  the 
quantum  efficiency.  This  can  be  calculated  approximately  by  means 
of  the  formula  below. 

Emu  -  [e.(TJ  +  EW(T„)  -  E,(-|  T)  -  EV(J>  T)]  X  0.409  (38) 

IV.  Numerical  calculations 


For  the  pseudo-one-dimensional  nonequilibrium  flow  equations, 
we  followed  the  method  used  in  [7] ,  and  adopted  the  difference  scheme 
given  in  [21] .  It  is  fairly  easy  to  show  that  this  computational 
scheme  possesses  second-degree  accuracy.  The  amount  of  computation 
is  much  smaller  than  that  required  when  a  Taylor  series  expansion 
with  second  degree  terras  is  used.  The  method  also  lends  itself  well 
to  computer  programming. 


A  problem  that  exists  in  computing  for  the  pseudo-one-dimensional 
flow  using  a  time-dependent  method  is  that  after  being  brought  to 
the  steady  state,  p  and  u  may  not  necessarily  satisfy  the  condition 
of  mass  conservation.  Originally,  from  equation  (25) ,  we  have 

dp  1_  djpmA') 

di  A 

After  the  steady  state  has  been  reached,  ^ 

puA  is  a  constant  at  all  points  of  the  nozzle.  However,  problems 
inherent  in  the  numerical  computation  often  cause  the  mass  flux  of 
the  convergent  solution  to  be  different  at  different  positions  x. 


dx 


(39) 

,  which  means  that 


12 


v.y 


vrv. 
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According  to  [7],  the  difference  in  the  values  of  puA  at  different 
portions  of  the  nozzle  can  be  as  high  as  20-30%  .  We  know  from  equa¬ 
tion  (39)  that 

Thus,  the  difference  in  the  mass  fluxes  at  neighboring  mesh  points 
is  not  only  dependent  on  the  difference  ratio  Ap/At,  but  also  on 
the  space  step  size  Ax.  To  reduce  the  size  of  A(puA),  besides 
greatly  increasing  the  number  of  computations  so  as  to  absolutely 
reduce  J Ap/At |,  one  may  resort  to  the  possibly  more  effective  way  of 
increasing  the  density  of  the  mesh  points.  Furthermore,  substituting 
the  right  hand  side  of  equation  (39)  for  the  expanded  right  hand 
side  of  equation  (25)  will  also  eliminate  errors  due  to  computations 
of  the  difference  ratio,  and  make  A(puA)  smaller. 

We  ran  the  program  that  we  compiled  on  the  TQ-16  computer  made 
in  this  country.  To  examine  the  various  factors  that  affect  the 
accuracy  of  computation,  and  to  arrive  at  a  way  of  reducing  comput¬ 
ing  time  under  given  requirement  on  accuracy,  we  first  did  a  large 
amount  of  experimentation  with  numbers . 

From  1/3  to  1/2  of  the  time  of  computing  for  each  time  step  is 
spent  on  the  computation  of  the  relaxation  time.  The  relaxation 
time  is  mainly  dependent  on  T  which  does  not  vary  much  between  two 
neighboring  time  steps,  the  amount  of  variation  decreasing  with 
increasing  convergence.  Therefore,  we  were  able  to  appreciably 
reduce  computing  time  by  means  of  skipping  several  time  steps  before 
making  a  computation  for  the  relaxation  time  x. 

To  reduce  computing  time,  Anderson  carried  out  the  computation 
for  the  entire  nozzle  in  two  sections,  making  use  of  mesh  points 
with  two  different  space  step  sizes.  We  used  tri-sectional  computa¬ 
tion  with  variable  intervals,  and  carried  out  the  computation  in 
each  section  only  down  to  a  few  points  downstream  the  boundary  point. 
This  speeded  up  the  computation.  However,  there  is  some  discrepancy 
in  the  results  computed  for  AuA  at  the  boundary  points  from  the 


different  sections.  These  kinks  can  be  smoothed  out  by  means  of 

recalculations  using  the  converged  solutions  as  initial  values. 

Our  computations  show  that  such  smoothing  has  but  little  effect  on 

the  gain,  E  and  flow  parameters  downstream  from  the  nozzle  exit  that 
max 

we  are  mainly  interested  in. 

Computing  time  can  be  significantly  reduced  by  using  the  con¬ 
verged  solutions  for  a  nonequilibrium  flow  as  the  initial  distribu¬ 
tion  for  another  set  of  solutions  for  a  problem  with  slightly  diff¬ 
erent  conditions.  This  method  is  especially  useful  when  one  wishes 
to  examine  the  effects  of  continual  changes  in  a  certain  parameter 
(such  as  Tq  or  pQ,  etc.) . 

We  have  carefully  examined  the  effects  on  convergence,  accuracy 
of  the  results  and  computing  speed  of  the  choice  of  the  position  of 
the  boundary  point  between  different  sections,  the  position  of  the 
downstream  end  of  each  section  of  computation,  the  initial  value 
distribution,  the  density  of  mesh  points  in  every  section  of  the 
nozzle,  the  number  of  steps  one  skips  before  making  a  calculation  of 
x,  condition  for  convergence,  time  step  size,  and  the  position  of 
the  starting  point  of  the  subsonic  section  of  the  nozzle  (the  upstream 
end  of  the  entire  calculation) . 

It  is  worth  pointing  out  that  when  the  mesh  points  are  not  pro¬ 
perly  arranged,  even  though  the  non-steady-state  calculations  for 
the  nozzle  flow  can  still  converge,  there  often  appears  an  unsmooth¬ 
ness  in  the  subsonic  section.  The  parameters  have  fluctuating  values 
instead  of  values  that  vary  monotonically .  These  values  may  even 
become  greater  than  the  stagnation  values  (such  as  T,  p) .  Such  a 
phenomenon,  which  is  irrational  from  the  standpoint  of  physics,  is 
due  to  improper  numerical  procedures.  This  problem  can  usually  be 
overcome,  and  a  smooth  solution  can  be  obtained  by  increasing  the 
density  of  mesh  points,  especially  those  near  the  throat.  Of  course, 
this  lengthens  the  computing  time.  Luckily,  even  if  the  solution  is 
not  smooth,  the  resulting  parameters  downstream  from  the  nozzle  exit  are 
not  much  different  from  those  of  the  smooth  solution. 


In  addition,  improper  choices  of  the  time  step  size  will  also 
result  in  an  unsmooth,  or  even  divergent,  solution. 

V .  Results  and  discussion 

1.  Comparison  with  experimental  results.  To  examine  the 
effects  of  improvements  made  in  the  relaxation  model,  the  relaxation 
equations  and  the  data  used,  as  opposed  to  Anderson's  results,  we 
carried  out  computations  for  some  published  experimental  results. 

As  the  detailed  nozzle  shape-line  data  were  not  supplied  in  the 
references,  we  used  in  our  calculations  a  nozzle  of  the  smallest 
length  that  had  the  same  throat  height  and  approximately  the  same 
area  ratio.  Our  calculations  show  that  the  small  signal  gain  of  the 
nozzle  exit  is  not  sensitive  to  differences  in  the  shape-lines. 

(See  the  results  given  below) . 


To*1800°K,  pq=37 . 5  atm,  XCQ  =0.07,  Aexit/A*=50»  h*=0. 0356cm 

nozzle  of  minimum  length  O  experimental  data  [14] 

1 — computed  results  given  in  [13],  in  the  absence  of  Sharma  potential 
barrier;  2 — computed  results  given  in  [13],  in  the  absence  of  Sharma 
potential  barrier;  3 — computed  results  given  in  [7]  and  [8]  (old 
relaxation  rates);  4 — computed  results  given  in  [17]  (new  relaxation 
rates) ;  5 — computed  results  given  in  this  paper 


Figure  3.  Go~XH  0  curves  (l*27*cm  downstream  from  the  nozzle  exit) 

(A  comparison  of  difflrent  relaxation  models  and  relaxation  rates) 
To=1800°K,  P0=20  atm,  XC()2=0.07,  Aexit/A*=20'  h*=0.1 

nozzle  of  minimum  length  A  experimental  data  [14] 

1 — computed  results  given  in  this  paper;  2 — results  obtained  by 
using  the  relaxation  rates  given  in  [12]  and  the  relaxation  model 
given  in  this  paper;  3 — results  obtained  by  using  the  relaxation 
rates  given  in  [7]  and  [8]  (old  rates)  and  the  model  given  in  this 
paper;  4 — relaxation  rates  given  by  [17]  (new  rates)  and  the  model 
given  in  this  paper;  5 — relaxation  rates  given  in  [7]  and  [8]  (old 
rates)  and  Anderson's  model  (taken  from  [14]);  6 — relaxation  rates 
given  in  [17]  (new  rates),  and  Anderson's  model  (taken  from  [14]). 


It  can  be  seen  from  Figures  2  and  3  that  Anderson's  results  do 
not  agree  very  well  with  the  experimental  results  for  small  signal 
gains.  When  his  new  data  [17]  are  used,  the  results  computed  for 
the  case  with  A/A*=50  agree  fairly  well  with  experimental  results. 
The  results  for  the  nozzle  with  A/A*=20,  however,  are  much  lower 
than  the  experimental  results.  This  is  inexplicable  [14],  In  con¬ 
trast,  the  tendency  of  variation  of  our  computed  results  agree  with 
the  experimental  data.  The  theoretical  results  for  the  nozzle  with 
A/A*=20  agree  with  the  experimental  results,  while  those  for  the 
nozzle  with  A/A*=50  are  higher  than  the  experimental  results  by 
about  30%.  This  is  understandable  because  for  a  nozzle  with  greater 
expansions,  the  effects  of  neglecting  the  non-ideal  characteristics 
of  the  flow  and  the  two-dimensional  effects  in  the  computations,  and 
the  nonequilibrium  effects  of  vibration  in  the  design  of  the  nozzle 
are  aggravated.  The  effects  of  variations  and  deviations  in  the 
assembly  and  processing  of  the  nozzles  will  also  be  greater  (this  is 
almost  inevitable) .  Thus,  we  may  say  that  our  computations  have 


Figure  4.  Schematic  diagram  of  the  nozzle  shape-line 
(a)  nozzle  with  circular  throat;  (b)  nozzle  with  minimal  length 
1 — oncoming  stream;  2 — throat 

done  away  with  the  inconsistency  in  Anderson's  computations. 

It  can  be  seen  from  Figure  2  that  the  conclusion  that  Anderson 
made  from  his  procedure  of  computation,  which  was  based  on  old  data 
and  which  contained  relatively  large  errors,  that  increasing  the 
H20  content  of  second-generation  gasdynamic  lasers  will  keep  the 
gain  basically  the  same  is  not  correct.  Actually,  the  optimal  water 
content  of  gasdynamic  lasers  with  high  area  ratios,  high  stagnation 
temperatures  and  high  stagnation  pressures  is  still  around  1% ,  the 
gain  decreasing  rapidly  with  increasing  water  content. 

2.  Effects  of  the  nozzle  shape-line.  Supersonic  section:  We 
carried  out  calculations  for  a  nozzle  with  circular  throat  and  a 
nozzle  with  minimal  length  (Figure  4),  with  A/A*  ^  30  (M=5.0),  under 
similar  conditions.  The  results  are  given  in  Figures  5(a)  and  5(b). 
These  show  that  the  gain  of  the  nozzle  with  minimal  length  tends  to 
be  greater  by  less  than  8%  (even  less  for  Tn  >  1500 °K) ,  and  its  E__v 
is  greater  by  around  10% .  In  view  of  considerations  related  to  pro¬ 
cessing,  assembly  and  flow  fields,  it  would  be  all  right  to  give  up 
this  advantage  and  use  the  nozzle  with  circular  throat  instead. 

Subsonic  section:  In  the  calculations,  the  subsonic  section 
of  the  nozzle  is  made  up  of  a  circular  arc  and  a  straight  line 
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Figure  5.  Relation  between  (Emax)  exit,  <G0>exit  and  T(XC02) 

(CO  gauged  combustion)  po=30  atm,  xh2o=0*01 

- nozzle  of  minimum  length  - nozzle  with  circular  throat  shape - 

line 

tangent  to  this  arc  and  at  a  definite  angle  with  respect  to  the  cen¬ 
ter  line.  When  the  radius  of  the  arc  stays  the  same  while  the 
angle  of  inclination  a  varies  from  30°  to  75°,  both  Gg  and  Emax 
remain  basically  the  same.  The  Gn  at  the  nozzle  exit  varies  from 

0.796  to  0.799  (m  ),  and  E  changes  from  14.90  to  14.93  (J/g) . 

max 

When  the  angle  of  inclination  is  fixed,  Gn  and  E _ _  vary  with  R. 

u  max  l 

(Figure  6) .  It  can  be  seen  from  Figure  6  that  with  decreasing  radius 

of  the  arc  of  the  subsonic  section  near  the  throat,  both  G_  and  E 

u  max 

increase  monotonies lly,  and  reach  their  maximum  values  when  the  arc 
becomes  a  straight  line  (R^*0) .  The  range  of  variation  is  a  few 
percent.  Our  computations  show  that  the  part  of  the  subsonic  section 
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Figure  6 .  Relation  between 
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Figure  7.  Relation  between 

throat  height  and  Gn,  E  „ 

0  max 


for  which  the  area  ratio  is  greater 
than  5  there  is  almost  no  freezing 
effect,  and  its  shape-line  may  be 
chosen. 


Optical  cavity  section:  As 

the  vibrational  energy  stored  in  N2  and  C02  will  eventually  be 
passed  to  the  average  motion  via  relaxation  processes,  thus  raising 
the  temperature  of  the  average  motion  downstream,  this  effect  will 
be  more  pronounced  in  the  presence  of  laser  radiation.  To  increase 
the  gain  downstream,  from  the  optical  cavity  section,  one  may  allow  the 
optical  cavity  section  to  expand  further. 

3.  The  effect  of  the  nozzle  expansion  area  ratio.  Since  the 
relaxation  processes  transfer  the  vibrational  energy  to  the  average 
motion,  the  temperature  of  the  average  motion  of  the  nonequilibrium 
flow  will  be  higher  by  several  tens  of  degrees  than  that  computed 
for  an  isentropic  flow.  This  is  disadvantageous  for  bringing  about 
population  inversion.  In  order  to  fully  utilize  the  stored  vibra¬ 
tional  energy  to  increase  the  freezing  effect  of  the  nozzle,  one  may 
increase  the  area  ratio  to  reduce  the  temperature  downstream.  This 
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TABLE  2.  Relation  between  stagnation  pressure  and  G  (m  ) ,  E  (J/g) 

u  in  3.x 
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1 — pressure;  2 — G_  at  the  exit;  3--location  of  maximum  of  GQ; 

4 — maximum  G_;  5--E  at  the  exit;  6 — downstream;  7 — exit 
u  max 

will  result  in  a  significant  increase  in  both  G.  and  E  ,  and  is  a 

3  0  max' 

very  important  means  for  improving  gasdynamic  laser  performance. 

The  computed  results  are  given  in  Figures  5 (a)  and  5 (b) .  GQ  of 
the  nozzle  with  A/A*=50  is  30-40%  higher  than  that  of  the  nozzle 
with  A/A*=30,  while  Emax  of  the  former  is  5%  higher  than  that  of  the 
latter.  For  combustion  type  gasdynamic  lasers,  the  optimal  value  of 
the  stagnation  temperature  increases  slightly  with  increasing  area 
ratio . 

4.  Changes  in  the  throat  height.  Dimensionless  mesh  parameters 
are  used  in  nozzle  designs.  A  reduction  in  the  throat  height  is 
equivalent  to  a  shortening  of  the  blades.  This  means  that  a  shorter 
time  is  required  for  reaching  the  same  degree  of  expansion,  which  is 
obviously  advantageous  for  freezing  the  vibrational  energy  and  bring¬ 
ing  about  population  inversion  (Figure  7) .  Reducing  throat  height 

is  also  one  way  of  improving  the  performance  of  the  unit.  However, 
reduced  throat  height  and  blade  length  means  greater  difficulty  in 
blade  processing  and  assembly.  In  addition,  the  effect  of  the  bound¬ 
ary  layer  may  become  too  large  to  be  negligible. 

5.  Changes  in  the  stagnation  temperature.  In  combustion  type 
gasdynamic  lasers,  the  energy  source  for  the  thermal  excitation  comes 


from  the  combustion  of  the  fuel  (e.g.,  COg ,  I^) .  Therefore,  changes 
in  the  stagnation  temperature  can  be  brought  about  by  changes  in  the 
composition  of  the  gas  mixture.  No  computation  for  such  conditions 
with  practical  significance  has  been  found  in  the  literature.  We 
present  our  results  in  Figure  5. 

Generally,  GQ  has  a  maximum  near  1400-1600 °K  (the  value  being 
slightly  different  for  different  shape-lines) .  Emax  increases  in 
general  with  Tg,  but  the  tendency  is  not  very  strong  at  high  temper¬ 
atures.  Thus,  overly  increasing  the  operating  temperature  for  a 
given  fuel  will  not  be  very  helpful.  Furthermore,  higher  TQ  means 
greater  heat  loss  as  well  as  greater  difficulty  in  the  manufacture 
of  the  unit. 

6.  Changes  in  the  stagnation  pressure.  An  increase  in  pressure 

results  in  an  increase  in  the  density  of  particles,  more  frequent 

collisions  and  increased  relaxation  rate.  This  will  cause  GQ  and, 

especially,  EQ,  to  drop.  The  computational  results  are  as  shown  in 

Table  2  and  Figure  8 .  The  maxima  of  E  and  Gn  decrease  monotonic- 

max  0 

ally  with  increasing  pg ,  and  the  optimal  point  of  the  gain  moves  in 
the  upstream  direction.  All  these  are  disadvantageous  for  the  per¬ 
formance  of  the  unit. 

The  light  output  of  the  operating  medium  per  unit  mass  can  be 
increased  by  improving  the  starting  performance  of  the  unit  and  the 
performance  of  the  pressure  diffuser,  and  reducing  the  operating 
pressure.  In  contrast,  increasing  the  stagnation  pressure  will 
reduce  the  light  output  per  unit  mass,  but  will  increase  the  totalut 
output  power  of  a  unit  of  the  same  size. 

7.  Changes  in  the  composition.  The  effect  on  Tg  of  a  change 
in  the  composition  has  been  given  above.  In  combustion  type  gas- 
dynamic  lasers,  inevitable  heat  losses  result  in  the  consumption  of 
more  heat  than  expected  to  reach  the  desired  Tg,  thus  causing  Xc 

(or  X„  )  to  increase  slightly.  To  examine  this  effect,  we  performed 


TABLE  3.  Effect  of  composition  on  C- 
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Figure  8.  Effect  of  pressure  Figure  9.  Effect  of  X  _  on 
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the  calculations  shown  in  Table  3.  It  can  be  seen  from  this  table 
that  Gq  is  slightly  higher  for  the  gas  mixture  with  a  higher  com¬ 
position  of  the  fuel.  This  is  because  Gq  is  directly  proportional 
to  Xc  (see  Equation  (37)).  However,  an  increase  in  Xc  will  speed  up 
the  release  of  the  vibrational  energy,  reduce  the  effective  storage 
of  the  vibrational  energy  (XM  decreases) ,  and,  therefore,  cause  E 
to  drop  significantly.  Thus,,  the  amount  of  energy  stored  in  the 
molecules  that  can  be  used  for  lasing  is  reduced,  and  the  light  out¬ 
put  cannot  be  enhanced.  In  addition,  in  the  presence  of  laser  radia 
tion,  a  high  Xc  has  an  adverse  effect  on  the  saturation  effect  of 
the  gain.  Therefore,  in  examining  the  various  factors  affecting  the 
performance  of  the  unit,  one  should  not  just  consider  the  effect  on 
Gft  but  should  make  an  overall  analysis. 


To  obtain  the  same  stagnation  temperature,  one  can  change  the 
composition  of  the  gas  mixture.  As  the  heat  of  combustion  and  the 
specific  heat  of  the  products  are  similar  for  CO  and  one  may 
regard  gas  mixtures  with  constant  Xc+XH  as  having  approximately  the 
same  stagnation  temperature.  The  results  of  our  computation  are 
given  in  Figure  9 .  The  optimal  water  content  of  combustion  type  gas- 
dynamic  lasers  is  also  in  the  neighborhood  of  1% .  An  increase  in  the 
water  content  will  cause  G.  and  E  to  drop  sharply.  This  is  because 
increased  water  content  will  greatly  increase  the  radiationless  dis¬ 
sipation  of  vibrational  energy. 

8.  Effects  of  relaxation  model,  equations  and  relaxation  data. 

In  [7],  [17]  and  [12],  the  logarithm  of  the  relaxation  time  is  con- 

-1/3 

sidered  to  have  a  linear  relation  with  T  '  .  However,  this  linear 

relation  does  not  hold  very  well  in  the  greater  temperature  range 
that  includes  the  lower  temperatures.  Hence,  the  use  of  linear 
fitting  may  solve  some  of  the  problems  while  missing  out  on  others. 

In  this  paper,  we  have  taken  up  to  the  second  or  third  power  terms 
in  the  relation  between  the  relaxation  data  and  T-1,/3.  Therefore, 
strictly  speaking,  our  paper  is  different  from  the  above  mentioned 
papers  in  the  data  used  for  every  reaction.  We  list  only  the  major 
differences  in  Table  4.  It  can  be  seen  from  this  table  that  the  new 
rates  used  in  [17]  represent  changes  in  the  data  for  two  energy 
exchange  reactions  only,  and  the  fitting  to  experimental  the  relaxa¬ 
tion  data  is  still  not  very  good.  The  results  given  in  [12]  are 
closer  to  those  given  in  this  paper,  but  there  is  still  a  significant 
difference  in  the  rate  of  radiationless  dissipation  of  the  lower 
energy  level  of  the  laser. 

Figure  3  gives  the  results  obtained  (curves  1-4)  by  using  our 
relaxation  model  and  relaxation  equations,  and  several  different  sets 
of  relaxation  data.  (As  reaction  (3)  was  not  taken  into  consideration 
in  [7],  [17]  and  [12],  our  data  was  used  for  this  reaction).  It 
also  gives  the  results  taken  from  [1]  which  were  obtained  by  using 
the  model  and  equations  as  well  as  the  new  and  old  relaxation  rates 


given  in  [17]  and  [7]  (curves  5  and  6) .  Obviously,  there  exist  a 
large  difference  between  the  results  obtained  by  using  different 
sets  of  relaxation  data  with  the  same  relaxation  model  and  equations 
This  point  is  also  borne  out  by  the  computation  given  in  [13] .  In 
our  opinion,  even  the  new  data  given  in  [17]  are  not  adequate  enough. 
A  comparison  of  curves  3  and  5,  and  curves  4  and  6  shows  that  a  large 
difference  also  results  from  using  the  same  relaxation  data  but  diff¬ 
erent  models  and  equations.  The  resulting  changes  may  be  very  irreg¬ 
ular.  We  may  thus  deduce  that  the  errors  brought  about  by  the  relax¬ 
ation  model  and  equations  used  in  [7]  and  [17]  were  for  the  most  part 
offset  by  those  due  to  the  relaxation  data.  However,  as  the  cancel¬ 
lation  cannot  always  be  complete,  one  has  the  irrational  phenomena 
of  curve  6  being  much  lower  than  the  experimental  result. 

More  refined  examination  and  analysis  may  be  made  of  the  effects 
of  relaxation  model,  relaxation  equations  and  relaxation  data  on 
computational  results  for  nonequilibrium  flows. 


vavui.'.  /. 
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